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ABSTRACT 

We examine a large set of synthetic quasar spectra to realistically assess the potential 
of using the relative sizes of highly ionized near-zones in the Lya and Ly/3 forest as a 
probe of the neutral hydrogen content of the intergalactic medium (IGM) at z > 6. The 
scatter in the relative near-zone size distribution, induced by underlying fluctuations 
in the baryonic density field and the filtering of ionizing radiation, is considerable 
even for fixed assumptions about the IGM neutral fraction. As a consequence, the 
current observational data cannot distinguish between an IGM which is significantly 
neutral or highly ionized just above z = 6. Under standard assumptions for quasar ages 
and ionizing luminosities, a future sample of several tens of high resolution Lya and 
Ly/3 near-zone spectra should be capable of distinguishing between a volume weighted 
neutral hydrogen fraction in the IGM which is greater or less than 10 per cent. 
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1 INTRODUCTION 

Quasar absorption spectra are currently the premier obser- 
vational probe of the hydrogen reionization epoch. The aver- 
age amount of Lyman series absorption observed in quasar 
spectra is consistent with an intergalactic medium (IGM) 
with a volume weighted neut ral hydrogen fraction in excess 
of (.f m)v ~ 10" 3 ' 5 at z ~ 6 (|Fan et al.ll2006bl ; iBecker et all 
2006). However, obtaining more stringent constraints on the 
IGM neutral hydrogen fraction using this technique at z > 6 
is extremely difficult; larger neutral fractions produce sat- 
urated Lyman series absorption over substantial regions in 
the spectra of t hese quasars, the observa tional manifestation 
of which is the lGunn fc Peterson! (|l965l ) trough. 

Consequently, several techniques which may b e sensitive 
to larg er values of (/hi)v have been proposed (see lFan et al.l 
2006a and references therein). The sizes of transparent re- 
gions observed immediately blueward of the Lya emission 
line of z > 6 quasars, if equivalent to the size of the 
H ii regions surrounding the quasars, should be proportional 
to the ambient IGM neutral hydrogen fraction when adopt- 
ing assumptions for the quasar age and ionizing luminosit y 
(|Cen fc Haimanll2000l ; iMadau fc ReesfeoOOl ; lYu fc Lulliool) . 
Under these assumptions, these regi ons are consistent with 
(/hi)v > 0.1, jus t above z = 6 (jWvithe fc Loebl |2004 
IWvithe et al.ll2005l h although a recent independent analy- 
sis including modelling of the Gunn-Peterson trough damp- 



ing wing by iMesinger fc Haimanl (|2006T ) using many sim- 
ulated quasar sight-lines favours a slightly lower limit of 
</m)v > 0.033. 



However, recent numerical studies which correctly 
model the radiative transfer of ionizing photons around these 
quasars indicate that interpreting the sizes of these highly 
ionized Lya near- zones with respect to the neutral hydro- 
gen content of the IGM is complicated IIBolton fc Haehnelq 
12007 j ; iMaselli et aT]|2007l ; iLidz et all 1 20071 ). There are two 
main reasons for this, in addition to the uncertain quasar 
age and ionizing luminosity. Firstly, for small neutral frac- 
tions the Lya near-zone resembles the classical proximity 
zone of a luminous quas ar embedded in a highly ionized 
IGM jBaitlik et alJll988h and does not correspond to the 
extent of an H n region expanding into a substantially neu- 
tral IGM. The observationally identified Lya near-zone sizes 
can then substantially underestimate the size of the region 
of enhanced ionization surrounding the quasar, leading to 
an overestimate of the IGM neutral hydrogen fraction. Sec- 
ondly, even for fixed assumptions about the ionization state 
of the IGM, variations in the IGM density and the ioniz- 
ing background along different quasar sight-lines combined 
with radiative transfer effects lead to a considerable scat- 
ter in the near-zone sizes. Drawing robust constraints on 
(/hi) v from only a small s ample of spectra is therefore dif- 
ficult. iBolton fc Haehneltl l|2007al ) found the observed near- 
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zone sizes at z > 6 are consistent with neutral hydrogen 
fractions as small as {/hi)v = 10 -3 ' 5 . 

Nevertheless, further progress may be possible by also 
considering the size of the corresponding Ly/3 near-zone. 
The edge of the Ly/3 near-zone should trace the position 
of the quasar H n ionization front (IF) to smaller volume 
averaged neutral hydrogen fractions in the ambient IGM 
relative to the Lya near-zone. The difference in the sizes 
of these regions thus holds add itional information on the 
IGM neutral hydrogen fra ction dMesinger fc Haimanll20o3 ; 
iBolton fc Haehneltl l2007al . hereafter BH07a). In principle, 
the ratio of Ly/3 to Lya near- zone sizes should increase from 
unity to a maximum value of ~ 2.5 as the IGM neutral 
hydrogen fraction decreases (BH07a). In practise, however, 
the ratio can be much smaller than 2.5 even for very small 
neutral fractions and it will exhibit a large scatter due to 
fluctuations in the IGM density along the line of sight. 

In this letter we extend the work of BH07a to examine 
the use of Lya and Ly/3 near-zones as a probe of {/hi)v- 
Our work closely follows on from the concepts discussed 
in BH07a, although the analysis presented here differs in 
three important ways. We analyse several hundred differ- 
ent simulations of the radiative transfer of ionizing photons 
around z > 6 quasars, greatly improving the synthetic spec- 
tra statistics. An approximate treatment of an inhomoge- 
neous ionizing backgroun d, expected towards the tail end 
of hy drogen reionization (|Fan et al.ll2006bl ; IWvithe fc Loebl 
2006) is now also included within the radiative transfer sim- 
ulations. Lastly we adopt a more robust method for defining 
the sizes of both Lya and Ly/3 near- zones. As a consequence, 
we are able to more realistically assess the potential of using 
quasar near-zones as a probe of the IGM neutral hydrogen 
fraction beyond z = 6. 



2 SIMULATIONS 

2.1 Radiative transfer implementation and the 
IGM density distribution 

We use the one dimensional, multi-frequency photon- 
conserving algorithm described and tested in BH07a to 
compute the radiative transfer of ionizing photons around 
quasars. This is combined with density distributions drawn 
from the 400 3 hydrodynamical simulation of BH07a, ru n 
using the parallel TREE-SPH code GADGET-2 (ISpringelll2005l) . 
The ten most massive haloes in the simulation at z = 6.25 
were identified and baryon density distributions were ex- 
tracted in different orientations around them. Continuous 
density distributions 150/i _1 comoving Mpc in length were 
then constructed using the halo density distributions com- 
bined with other sight-lines drawn randomly from the sim- 
ulation volume. One end of the density distributions always 
lies at the centre of one of the identified haloes, where the 
quasar is assumed to reside during the radiative transfer cal- 
culation. This process was repeated to produce 240 unique 
sight-lines. 

Synthetic Lya and Ly/3 spectra are constructed from 
the output of the radiative transfer simulations following 
standard procedure (e.g. BH07a). The raw synthetic spectra 
are processed to resemble data obtained with the Keck tele- 
scope Echelle Spectrograph and Imager (ESI). The spectra 



are convolved with a Gaussian with a full width half max- 
imum (FWHM) of 66 km s _1 and rebinned onto pixels of 
width 3.5 A (R ~ 2500). Gaussian distributed noise is then 
added with a total signal-to-noise ratio of 20 per pixel at the 
continuum level and a constant read out signal-to-noise of 
80 per pixel. We also process some of the spectra to resemble 
data taken with the Keck telescope High Resolution Echelle 
Spectrometer (HIRES). These are convolved with a Gaus- 
sian with a FWHM of 7 km s _1 and rebinned onto pixels of 
width 0.25 A (R ~ 35000). Noise is then added with a total 
signal-to-noise ratio of 20 per pixel at the continuum level 
and a constant read out signal-to-noise of 100 per pixel. 

Random sight-lines were also drawn from the simulation 
at z = 5.12 to model the foreground Lya forest superim- 
posed on the Ly/3 absorption at z — 6.25. These were spliced 
together to produce a z = 5.12 Lya spectrum with the 
same wavelength coverage as the Ly/3 spectrum at z = 6.25. 
Lastly, the foreground Lya optical depths are rescaled (e.g. 
IBolton fc H aehnelt 2007b]) to reproduce the obser ved mean 
flux of the Lya forest at z = 5.12, (F) = 0.149 (|Songailal 
|2004) . 



2.2 Initial conditions and the fluctuating ionizing 
background model 

We run a total of 480 radiative transfer simulations arranged 
into two groups of 240 simulations. The density distributions 
for both groups have initial neutral hydrogen fractions dis- 
tributed uniformly in the range — 4< log(/Hi)v < 0, set by 
assuming the IGM is in ionization equilibrium with an ion- 
izing background with a power-law spectral index, Ob = 3, 
below the Lyman limit. For the first group of 240 simulations 
the ionizing background is assumed to be spatially uniform, 
as in BH07a. 

However, a spatially uniform ionizing background is 
likely to be a poor approximation at the tail end of reion- 
ization and this assumption may have a n impact on the 
near-zone sizes derive d from simulations l|Lidz et al.ll2007l ; 
I Alvarez "fc~ Abel 2007J). Therefore, for the second group of 
simulations we construct spatially fluctuating ionizing back- 
grounds using observationally determined Lyman break 
galaxy (LBG) lumino sity functions at high redsh ift. We do 
this using the model of lBolton fc Haehneltl (|2007bh . to which 
we refer the reader for further details. An example of the 
fluctuating photo- ionization rate per hydrogen atom, Thi, 
along a single sight-line at z — 6.25 is shown in Figure [1] 
The background value of Thi due to galaxies is several or- 
ders of magnitude larger in overdense regions due to the 
clustering of ionizing sources, although it is still dominated 
by the quasar radiation field within 10 proper Mpc of the 
host halo. In order to reproduce different initial values of 
(/hi)v in our simulations we rescale the mean of the fluctu- 
ating ionization rate along each sight-line by the appropriate 
amount. Note, however, that this simple model onl y applies 
durin g the post-overlap phase of reionization (cf. iGnedinl 
2000) when the ionizing photon mean free path is larger than 
the typical separation between ionizing sources. For a mean 
free path which is shorter than the typical ionizing source 
separation, large values of (/m)v can also be produced by 
the presence of small, hig hly ionized regions surrounded by 
an entirely neutral IGM (|Lidz et al 1 12007I ; Iai varcz fc Abell 
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Figure 1. Top: Example of the spatially fluctuating hydrogen 
ionization rate along a simulated quasar sight-line. The quasar 
is situated at the right hand side of the diagram. The dotted 
line shows the ionizing background due to galaxies only and the 
solid line also includes the contribution from the quasar computed 
using the radiative transfer implementation. Bottom: The corre- 
sponding baryonic overdensity along the sight-line. The ionizing 
background is enhanced in overdense regions due to the clustering 
of ionizing sources. 
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Table 1. The Lya and Ly/3 near-zone sizes around the five 
quasars with published spectra at z > 6.1. The sizes have been 
measured from the spectra using the definition given in Sec- 
tion |3T] 



Quasar R a [proper Mpc] Hp [proper Mpc] Rp/R a 



J1509 


- 1749 


6.4+1.2 


12.7 + 1.2 


2.0 


J1250 


+ 3130 


13.6 + 1.2 


9.3 ± 1.2 


0.7 


J1623 


+ 3192 


2.6 + 1.0 


3.1 ± 1.0 


1.2 


J1030 


+ 0524 


4.8 + 1.2 


6.0 ± 1.2 


1.2 


J1148 


+ 5251 


6.2 + 1.4 


5.6 ± 1.4 


0.9 



tion where the last pixel in an unsmoothed spectrum drops 
below a normalised flux threshold of F = 0.1. Although 
providing a better measure of the Ly/3 near-zone size, this 
definition still breaks down once the IGM as a whole be- 
comes highly ionized and the edge of the near-zones become 
ambiguous. To overcome these difficulties, we add the extra 
condition that the last pixel at which the flux drops below 
F = 0.1 must be followed by a gap of Az > 0.1 where the 
normalised flux remains below the F = 0.1 threshold. Any 
spectra which do not meet this criterion are considered to 
have unidentifiable Lya and Ly/3 near-zones and are rejected 
from our analysis. This occurs for only one and three sight- 
lines in the synthetic ESI and HIRES samples, respectively. 



120071 ). We discuss the possible impact of these pre-ionized 
regions on near-zone sizes later in this letter. 

The fiducial quasar parameters adopted in the radia- 
tive transfer simulations are an ionizing photon production 
rate of N — 2 x 10 57 s _1 , an age of £q = 10 7 yrs and 
a power-law spectral index of a s — 1.5 below the Lyman 
limit. The value of N is consistent with the z > 6 quasar lu- 
minosities infe rred by Fan et al,l 1120060) a nd independently 
constrained bv lMesinger fc Haimanl (|2006l ). Quasar ages are 
rather more uncertain, although 10 7 y rs is in agree ment 
within current obse rvational constraints (|Martini||2004l ') and 
theoretical models (Hopk ins et af]|2006j ). The impact of dif- 
fering values of N and £q on the sizes of quasar near-zones 
is discussed in detail in BH07a. 



3 RESULTS 

3.1 A consistent definition for Lya and 
Ly/3 near-zone sizes 

There are considerable variations in Lya and Ly/3 near-zone 
sizes even for fixed values of (/hi)v due to underlying fluc- 
tuations in the IGM density field. Therefore, a consistent 
definition for the sizes of the Lya and Ly/3 near-zo nes is very 
import ant when comparing individual spectra. iFan et all 
l|2006bh smoothed the observed Lya near-zone spectra to 
a resolution of 20 A and identified the near- zone sizes at the 
first pixel which dropped below a normalised flux threshold 
of F = 0.1. However, this definition is not suitable for the 
Ly/3 near-zone, where the transmission is much weaker and 
more patchy due to the foreground Lya forest. Alternatively, 
BH07a measured Lya and Ly/3 near-zone sizes at the posi- 



3.2 The impact of ionizing background 

fluctuations on the sizes of near-zones 

The sizes of the Lya and Ly/3 near-zones measured from our 
480 synthetic ESI spectra are shown as the red diamonds in 
Figure [2] The top panels show the sizes measured assum- 
ing a uniform ionizing background, while the bottom panels 
correspond to the sizes measured with the assumption of 
a spatially fluctuating ionizing background. We find a sub- 
stantial scatter in the near-zone sizes even for fixed assump- 
tions for (/hi)v- The scatter is rather similar for the spectra 
computed with either a uniform or spatially fluctuating ion- 
izing background, suggesting that differences in the density 
distribution and radiative transfer effects play a more im- 
portant role in generating this variation. There is also little 
evidence that enhanced ionization in the overdense regions 
increases the average near-zone size. Note, however, that 
during the overlap-phase spatial fluctuations in the ionizing 
background are expected to be larger than we have assumed 
here and can extend over wider scales. This is due to the 
rapid increase of the mean free path for ionizing photons in 
the short phase when the ionized regions percolate. 

The filled black circles with error bars in Figure [2] cor- 
respond to the average near-zone sizes in bins of width 
A[log(/Hi)v] = 0.5. The error bars correspond to the la 
standard error of the mean. There are three main regimes 
in the Lya near-zone size as a function of log{/Hi)v for 
the quasar age and luminosity adopted here. The first is 
when R a oc (fm)^ 1 ^ 3 for (/hi)v > 0.1. In this regime the 
Lya near-zone edge closely corresponds to the H n IF. The 
second regime is when the near-zone corresponds to a classi- 
cal proximity zone and the size is independent of the neutral 
fraction; in this case the H 11 IF actually lies ahead of the 
Lya near-zone edge but the residual neutral hydrogen be- 
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Figure 2. Top: The red diamonds correspond to Lya (left panel) 
and Ly/3 (right panel) near-zone sizes measured from 240 syn- 
thetic ESI spectra using the method described in Section I3.ll 
The data are plotted as a function of the volume weighed neu- 
tral hydrogen fraction in the IGM. A spatially uniform ionizing 
background has been assumed. The filled black circles with ltr 
error bars correspond to the average near-zone sizes in bins of 
width A[log(/m)v] = 0-5- Bottom: As for the top panel except 
assuming a spatially fluctuating ionizing background. The filled 
blue and cyan squares in all panels correspond to near-zone sizes 
we derive from the published observational data, given in Table[T] 



hind it is large enough to produce saturated Lya absorption. 
Lastly, for (fm)v < 10 -3 ' 5 the measured Lya near-zone 
sizes increase again due to additional transmission originat- 
ing from regions now highly ionized by the ionizing back- 
ground. The Ly/3 near-zones exhibit a similar trend to the 
Lya data, except due to the smaller Ly/3 absorption cross- 
section they are able to trace the position of the H n IF to 
lower IGM neutral hydrogen fractions, reaching Rp oc (/hi)v 
at ( /hi)v — 10~ 2 . Therefore, as noted by BH07a (but see 
also Mesinger fc Hai man 2004 for a different interpretation) , 
evidence for Rp/R a > 1 may provide an interesting con- 
straint on (/hi)v- 

Finally, we plot the Lya and Ly/3 near-zone sizes we 
measured from the five quasars at z > 6.1 with published 



spectra as filled blue and cy an squares in Figure [2l (|Fan et al.l 
l2006bl ; IWillott et ai1l2007n . We have omitted J1030 + 0524 
at z — 6.2 from this sample since this is a broad absorption 
line quasar which c omplicates the me asurement of the near- 
zone size. Following iFan et al.l l|2006bl ). the sizes are rescaled 
to a common absolute magnitude of M1450 = —27 by as- 
suming the near-zone size is proportional to N 1 ^ 3 , although 
TV 1 / 2 may be more appropriate (BH07a). This magnitude 
corresponds to N ~ 1.9 x 10 57 s _1 for the quasar spectrum 
adopted in this work. The lower limit of (/hi)v > 10 5 
measured from the Gunn-Peterson trough limits is assumed. 
An uncertainty of Az — 0.02 in the systemic redshift of the 
quasar is adopted when determining the systematic error 
on the near-zone sizes. We do not show the systematic un- 
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Figure 3. Left: The ratio of Ly/3 to Lya near-zone sizes as a func- 
tion of the volume weighted neutral hydrogen fraction, measured 
from synthetic ESI spectra constructed assuming a spatially fluc- 
tuating ionizing background. The red diamonds show the individ- 
ual measurements and the filled black circles with la error bars 
correspond to the averages in bins of width A[log(/ni)v] = 0.5. 
The filled blue and cyan squares correspond to the near-zone sizes 
measured from the published data, given in Table [l] Right: As for 
the left panel, except the data are now measured from synthetic 
HIRES spectra. 



certainties on the measured near-zone sizes in Figure [2] for 
clarity, but the data are given in Table [1] The observational 
data exhibit a large scatter similar to that seen in the syn- 
thetic data, and are consistent with a wide range of neutral 
fractions. 



3.3 Probing the neutral fraction with the ratio of 
Ly/3 and Lya near-zone sizes 

A potentially interesting constraint on (/hi)v may be ob- 
tained by considering the ratio of the Ly/3 to Lya near- 
zone sizes. The red diamonds in the left panel of Figure [3] 
correspond to Rp/R a as a function of log(/m)v for the syn- 
thetic ESI spectra constructed using the fluctuating ionizing 
background model. The filled black circles show the average 
values of Rp/Ra in bins of width A[log(/Hi)v] = 0.5. The 
error bars correspond to the la standard error of the mean. 
The earlier analysis by BH07a found R@ / R a > 1 was consis- 
tent with (/hi)v ^ 10 -2 . However, this was based on a small 
sample of 20 synthetic spectra which used only four different 
initial values of (/hi)v- The improved statistics and wider 
sampling of (/hi)v in this work indicate Rp/R a > 1 when 
(/m)v^0.1 at the 3er level for a sample of around 30 ESI 
spectra. 

Assuming that Rp/R a > 1 is indeed consistent with 
(/hi)v^0.1, how reliable a constraint on (/hi)v can the 
current observational data provide? The o bserved near-zone 
sizes we measure fr om the published data (|Fan et al.ll2006b1 ; 
IWillott et aT]|2007l ) at z > 6.1 are plotted as the filled blue 
and cyan squares in Figure [3] Once again, the observational 
data are consistent with a wide range of (/hi)v- The mean 
of the observational data favours a value of Rp / R a — 1.2, al- 
though with only five data points this result is certainly not 
a significant one. The simulated data indicates that about 
20 ESI spectra would be required for a 2a result. Obtaining 
reliable constraints on (/hi)v using this technique with the 
current data is probably not yet possible. 

Higher resolution quasar spectra would enable a more 



© 0000 RAS, MNRAS 000, 000-000 



A closer look at quasar near-zones 5 



accurate determination of the sizes of the Lya and Ly/3 near- 
zones. This is demonstrated in the right panel of Figure[3]in 
which our simulated spectra have been processed to resemble 
data taken with HIRES. The scatter in the data is signif- 
icantly reduced with the higher resolution spectra. Strong 
evidence for a neutral fraction {/hi)v < 1CP 1 should be ob- 
tainable with a sample of 30 high-resolution spectra at the 
3 — 5cr level. 



3.4 The effect of quasar age, ionizing luminosity 
and surrounding galaxies 

An important caveat to the argument presented in this work 
is that the (/hi)v threshold where Rf)/R a > 1 becomes 
statistically significant is also degenerate with the age of 
a quasar, its ionizing luminosity and the size of any pre- 
existing ionized region around the quasar host halo. The 
observable Rp/R a > 1 corresponds to the (/hi)v thresh- 
old at which -Rhii > R a , where 7?hii is the extent of the 
H n IF. The position of the H n IF will scale approximately 
as i?Hii oc (fm)^ 1 ^ 3 N 1,/3 tQ 3 but the saturated Lya near- 
zone size will scale as R a oc N 1/2 only (BH07a). Therefore, 
if the mean quasar age is ten times larger than the fiducial 
10 7 yrs, the H n IF will be around twice as far from the 
quasar but the Lya near-zone size will remain unchanged. 
The (/hi)v threshold corresponding to Rp/ R a > 1 will thus 
be ten times larger. Alternatively, a smaller quasar age will 
lower the (/hi)v threshold corresponding to R/s/R a > 1. 
Disentangling this quasar age and neutral fraction degen- 
eracy from the data is potentially very difficult without a 
good independent estimate of the typical quasar age. Note, 
however, the degeneracy with N is less dramatic, since both 
the near-zone size and the H u IF position scale similarly 
with N. Finally, luminous quasars are expected to be hosted 
by rather massive dark matter haloes and the reionization 
process may be further advanced in their immediate envi- 
ronment compared to a typical region of the Universe. If 
pre-existing ionized regions produced by clustered galaxies 
lower the local value of (/hi)v around these quasars, the 
IFs will tr avel further than one expects for a unifo rmly ion- 
ized IGM (|Lidz et al.ll2007l ; lAlvarez fc Abelll2007h . The ob- 
served near-zone sizes may then exhibit Rp / R a > 1 even if 
(/hi)v > 0.1. 



4 CONCLUSIONS 

We have carefully examined the potential of the ratio 
of Ly/3 to Lya near-zone sizes as a probe of the IGM 
H i fraction at z > 6 with a large set of detailed radiative 
transfer simulations. Adopting a robust, simultaneous defini- 
tion for Lya and Ly/3 near-zone sizes, we find that a sample 
of around 30 ESI quasar spectra would be required to dis- 
tinguish between an IGM which has an H i fraction greater 
or less than 10 per cent at the 3cr level. Unfortunately, the 
current observational data, consisting of five quasar spectra 
at 2 > 6.1, are too few to place significant constraints on 
(/hi)v- Our results also suggest that, at least in the post- 
overlap phase of reionization, fluctuations in the ionizing 
background will have a negligible effect on near-zone sizes. 
However, if the typical quasar age is substantially longer 
than the canonical 10 7 yrs, or if pre-existing ionized regions 



surround quasar host haloes during the pre-overlap stage 
of reionization, it may still be very difficult to distinguish 
between a highly ionized or substantially neutral IGM. At- 
tempts to disentangle the former effect from the data may 
benefit from alternative statistics, such as the detection of a 
Gunn-Peter son trough damping wing l|Mesinger fc Haimanl 
|2004 | 2006), while distinguishing the latter effect may be 
possible with a sample of z > 6 gamma ray burst (GRB) 
spectra, which should probe less biased regions compared 
to luminous quasars. However, only one such spectrum has 
been obtained so far |Totani et al.|[2006ri and its interpreta- 
tion is hampered by damped Lya absorption from neutral 
hydrogen in the GRB host galaxy. 

High-resolution data significantly improves the ability 
to constrain the neutral fraction, and HIRES spectra o f 
2 > 6 quasars are already in existence dBecker et al.ll2005l ). 
Next generation infra-red surveys fsee lLawrenceH2007l for a 
recent review) such as the UKIRT Infrared Deep Sky Survey 
(UKIDSS) are expected to find s everal tens of new quasars 
at 2 ~ 6 jVenemans et al]|2007t ). Future data sets should 
thus be able to provide an interesting limit on the IGM neu- 
tral hydrogen fraction from the relative sizes of Lya and 
Ly/3 near-zones. 



REFERENCES 

Alvarez, M. A. fc Abel, T. 2007, MNRAS submitted, 
|astro-ph/0703740| 

Bajtlik, S., Duncan, R. C, & Ostriker, J. P. 1988, ApJ, 
327, 570 

Becker, G. P., Rauch, M., fc Sargent, W. L. W. 2006, ApJ 
submitted, astro-ph/0607633 

Becker, G. D., Sargent, W. L. W., & Rauch, M. 2005, in 
IAU Colloq. 199: Probing Galaxies through Quasar Ab- 
sorption Lines, Williams, P., Shu, C.-G. & Menard, B. 
eds., p. 357 

Bolton, J. S. & Haehnelt, M. G. 2007a, MNRAS, 374, 493 
Bolton, J. S. fc Hae hnelt, M. G. 2007b, MNRAS submitted, 
|astro-ph/0703306| 

Cen, R. & Haiman, Z. 2000, ApJL, 542, L75 
Fan, X., Carilli, C. L., & Keating, B. 2006a, ARA&A, 44, 
415 

Fan, X. et al., 2006b, AJ, 132, 117 

Gnedin, N. Y. 2000, ApJ, 535, 530 

Gunn, J. E. & Peterson, B. A. 1965, ApJ, 142, 1633 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 
Robertson, B., & Springel, V. 2006, ApJS, 163, 1 

Lawrence, A. 2007, A&G in press. larXiv:0704.0809l 

Lidz, A., McQuinn, M., Zaldarriaga, M., Hernquist, L., & 
Dutta, S. 2007, ApJ submitted, |astro-ph/0703667| 

Madau, P. & Rees, M. J. 2000, ApJL, 542, L69 

Martini, P. 2004, in Coevolution of Black Holes and Galax- 
ies, Ho, L. C. ed., p. 169 

Maselli, A., Gallerani, S., Ferrara, A., & Choudhury, T. R. 
2007, MNRAS, 376, L34 

Mesinger, A. & Haiman, Z. 2004, ApJL, 611, L69 

Mesinger, A. & Haiman, Z .2006, ApJ, 660, 923 

Songaila, A. 2004, AJ, 127, 2598 

Springel, V. 2005, MNRAS, 364, 1105 

Totani, T., Kawai, N., Kosugi, G., Aoki, K., Yamada, T., 
Iye, M., Ohta, K., & Hattori, T. 2006, PASJ, 58, 485 



© 0000 RAS, MNRAS 000, 000-000 



6 J.S. Bolton & M.G. Haehnelt 



Venemans, B. P., McMahon, R. G., Warren, S. J., 
Gonzalez-Solares, E. A., Hewett, P. C, Mortlock, D. J., 
Dye, S., & Sharp, R. G. 2007, MNRAS, 376, L76 
Willott, C. et al. 2007, AJ submitted, arXiv:0706:0914 
Wyithe, J. S. B. & Loeb, A. 2004, Nature, 432, 194 
Wyithe, J. S. B. & Loeb, A. 2006, ApJ, 646, 696 
Wyithe, J. S. B., Loeb, A., & Carilli, C. 2005, ApJ, 628, 
575 

Yu, Q. & Lu, Y. 2005, ApJ, 620, 31 



© 0000 RAS, MNRAS 000, 000-000 



